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Abstract 

The one— and two- particle densities of up to four interacting electrons with spin, 
confined within a quasi one-dimensional "quantum dot" are calculated by numerical 
diagonalization. The transition from a dense homogeneous charge distribution to a 
dilute localized Wigner-type electron arrangement is investigated. The influence of 
the long range part of the Coulomb interaction is studied. When the interaction is 
exponentially cut off the "crystallized" Wigner molecule is destroyed in favor of an 
inhomogeneous charge distribution similar to a charge density wave . 

PACS. 71.45.Lr - Charge-density-wave systems. 

PACS. 73.20.Dx - Electron states in low-dimensional structures. 

PACS. 73.20.Mf - Collective excitations. 
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Single charges dominate the electronic properties of submicron structures at low tem- 
peratures. Due to small capacitances the charging energy associated with the addition of 
one electron into a given structure can exceed the thermal energy. Two recently discovered 
important phenomena can be explained by single electron charging effects: the Coulomb 
blockade of the dc-current through small tunnel junctions fl], |2|, and the periodic oscilla- 
tions of the conductance of quantum dots [EL ffl . 

In contrast to metallic systems semiconductor nano-structures allow to reduce the num- 
ber of electrons in a quantum dot by changing a gate voltage. Optical absorption experiments 
were performed on quantum dots which contain only AT = 2 ... 4 electrons ||. For very low 
electron densities additional effects in the transport properties can be expected due to the 
increasing importance of the Coulomb interaction || since the electrons tend to "crystallize" 
into an inhomogeneous ground state J/]. Besides Coulomb effects, linear || and non- lin- 
ear || transport experiments show fine structure in the current voltage characteristics that 
can be traced back to the granularity of the charge density distribution and to the lowest 
collective excitations |10[, respectively. 

In this paper, we demonstrate that for a quasi-one dimensional (ID) system containing 
a few electrons the charge density q{x) becomes inhomogeneous for sufficiently large mean 
electron distance r s . Our results indicate that there are two consecutive regimes of charge 
localization. First, g(x) starts to become inhomogeneous with increasing r s . The correspond- 
ing energetically lowest excitations are associated with changes in symmetry and total spin 



lT|]. When r s is further increased, g(x) vanishes almost completely between well separated 
maxima. In analogy with the Wigner lattice, we denote this limiting configuration in a finite 
system as a "Wigner molecule". Two classes of elementary excitations can be identified 
in this limit: tunneling between different permutational arrangements of the electrons and 



vibrations of the charge density R ITT 



We consider N interacting electrons including spin, confined in a ID square well po- 
tential of finite length L. The height of the potential was assumed to be finite but large. 
The Hartree-Fock approximation is known |T2[ to overestimate the ferromagnetic state and 



cannot reproduce the true ground state which is antiferromagnetic in ID [|T^]. We calculate 
eigenvalues and eigenstates of the Hamiltonian exactly. 
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The kinetic energy 



^H(f) 2 E^tc n ^E H (f) 2 ^0 (1) 



and the Coulomb interaction 



Eh \J~J ^ Vnwsmn^i^cl^Cn^C^a = E B ( -^-J Hi (2) 



a 4 
a, a 1 



scale like 1/L 2 and 1/L , respectively. Thus the latter dominates at large L. The natural 
energy and length scales are given by the Hartree En = e 2 /a-Q , and the Bohr radius, 
a-Q = eh 2 /me 2 (e relative dielectric constant, m effective electron mass). c\ la creates a 
particle with energy e n (= (nn) 2 /2 for an infinitely high well) and spin a . V n4Tl3n2ni is 
the matrix element of the interaction potential V(x,x') = e 2 /eJ(x — x') 2 + A 2 . The cutoff 



at short distances simulates a small transversal spread A < I of the wave functions and 
leaves the V nin3n2ni finite. In the calculations we assumed X/L = 2 • 10~ 4 . The eigenvalues 
of H = Eu(aB I L)[(a B / L) H + Hi] depend only weakly on A when A 1 0. To investigate 
the effect of the long range part of the interaction, in some of the calculations an exponential 
cutoff at large distances was introduced. No charge neutralizing background was taken into 
account since it is not expected to modify the results for fixed iV qualitatively. 

The Hamiltonian matrix in the basis of the Slater determinants was diagonalized numer- 
ically. The latter were constructed from the single electron states f n {x)X(7 where ip n is a 
spatial function and Xc a spinor with a =|, j . Considering 1 < n < M leads to a matrix 

of the rank R = ( ^ ^ (R < 1.5 • 10 4 , for M = 10 ... 17 ). We used Lanczos procedures 

when R > 3 • 10 3 . 

The one- and two-particle densities were determined from the normalized eigenfunc- 
tions ip(xi, xjy) = Yjv=i (biz/v^/VT) det{(p nx (xj)xai}, were v denotes the configurations 
{rii<7i, ...,n7vO"jv} • The coefficients h v can be chosen real. In second quantization the one- 
particle density is given by 

^) = E(4 N) \^)^ a (x)\4 N) ) , (3) 

were \ipQ ) = Z^ = ib°|^\ is the iV-electron ground state. ^l(x) = J2n=i V 9 n( x ) c n. ( 7 creates 
an electron at position x with spin a and | z/) = c\ n(7i ■ ■ ■ 1 0) is a noninteracting Slater 
determinant. Straightforwardly we have 
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Q( x ) = h vK> H 9n{x)Vn>{x){v'\i >(T C n , a \v) . (4) 



v '=\ 



The development of maxima and minima in g{x) with increasing mean electron distance 
r s = L/(N — 1) is shown in figure 1 for iV = 3,4 . The finite height of the square well 
potential was assumed to be proportional to 1/L 2 such that the number of non-interacting 
bound states was independent of L. This implies a non-vanishing but small g(x) outside the 
box. 

We distinguish three characteristic regimes of electron densities. For small distances, 
r s <,0.1as , the Coulomb interaction perturbs Eu(o,b/ L) 2 H only weakly such that g(x) 
is dominated by the lowest occupied single particle states. This explains the minimum at 
x = (Fig.l). g(x) changes qualitatively at r s « a-Q- A structure consisting of N 
peaks emerges. The "critical" length for this transition is of the same order as found in 
] i] for the transformation from an almost non-interacting energy spectrum to the multiplet 



spectrum dominated by interaction. When r s increases further, say r s >, lOOae, g{x) vanishes 
almost completely in finite regions between the maxima indicating a fully established Wigner 
molecule. In this limit the ground state energy can be approximated by that of a chain of 
elementary charges at equal distances r s . 

In order to investigate the influence of the long range part of the Coulomb interaction, 



we introduced an exponential cutoff V(x,x') oc e~ a \ x ~ x '\ / y (x — x') 2 + A 2 , were 1/a is the 
cutoff distance. The inset of figure 1 shows g(x) for N = 3 and a = lOa^ 1 . Despite 
the extremely short range of the interaction pronounced maxima are observed. However, in 
contrast to a = the density between the maxima is much higher even for very large r s . 



The inset of figure 1 rather resembles a charge density wave JL4] . We conclude that the long 
range part of the Coulomb interaction is essential for the charge condensation and for long 
range density-density correlations. The ground state energy for a > r" 1 does not converge 
towards the Coulomb energy of iV elementary point charges for L — > oo . 

For large r s (r s ^> Ob) the spectrum consists of multiplets []15| |. Their energetic separation 
Q correspond to vibrational excitations. Their internal splitting A is due to tunnelling 
through the Coulomb barrier between adjacent permutational arrangements of the electrons 
[ [TT|j . Asymptotically, Q decreases as r~ 7 . Power law behavior is also obtained for a > . 
However, the exponent 7 changes from 7<,1.5 (a = 0) to 7 ~ 2 (a > 0) (figure 2). 
The vibrational excitations are independent of the symmetry of the spatial part of the N- 
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electron wave function and of the total spin. In analogy to the argument for a = leading 



to 1 < 7 < 1.5 [FT], the almost flat two-particle potential between adjacent charge density 
peaks for a ^> r" 1 leads to 7 = 2. Each electron is practically moving freely within a 
box of length r s . On the other hand, A decreases roughly exponentially with increasing L, 
even for a > (inset of figure 2). However, A is larger and the decrease is less rapid when 
r s is increased than for a = 0. The enhancement results from the reduced thickness of the 
screened Coulomb barrier such that the corresponding tunnelling integral (oc A) increases. 
The modified dependence on r s is due to the fact that only for r s < a^ 1 the thickness 
of the interaction barrier varies with r s . We find that the tunnelling excitations can be 
described in terms of correlated, localized basis functions even for a > . However, in order 
to obtain similar quantitative agreement with the numerically determined eigenvalues as for 
a = a larger mean electron distance r s is required to fulfill A -C . As a consequence 
the characteristic length for the transition from almost non-interacting to the interaction 
dominated behavior increases with increasing a . 

The pronounced correlations between the electrons in our model are also observed in the 
pair correlation function of the iV-electron state | ' 



Qc (x,X f ) = £ |tft {x) vfrt,^/) 9 ^ x) | V ,(iV) ) (5) 

g c is related to the density-density correlation function g 2 via Q2(x,x') = g(x)S(x — x') + 
g c (x,x') . Physically, g c is the probability to find an electron at position x when another 
electron sits at x' (jdxj dx' g c (x, x') = N(N — 1) ). Figure 3 shows a contour plot of g c for 
the N = 3 ground state (S = 1/2) with L = 9.5ae • Practically the same result is obtained 
for the spin polarized excited state (S = 3/2). The difference between £ c of the ground state 
and the spin polarized state is less then 1%, independent of x, and x' . 

It is known that the Hartree Fock approximation overestimates the tendency of electrons 
with parallel spins to avoid each other and underestimates the correlations between electrons 
with antiparallel spins. Therefore large differences between correlation functions for different 
spin configurations are expected in this approximation. Recently Pfannkuche et al. demon- 
strated this property explicitly for three electrons in a two dimensional parabolic potential 



|L6] . The weak sensitivity of the pair correlation function on S shows that highly correlated 
electrons in a "slim" quantum dot cannot be treated within molecular field approximation. 
In summary, we have analyzed one- and two-particle densities for the lowest eigenstates 
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of N = 2, 3, 4 electrons confined in a quasi-lD square well. Three different regimes were 
identified. At high electron densities the Coulomb interaction is only a weak perturbation. 
Between a^ 1 > r^ 1 >, lCT 2 ^ 1 the electrons tend to localize. Simultaneously the eigen- 
value spectrum develops a multiplet structure that is dominated by the interaction. Both, 
the inhomogeneity of the charge density and the multiplet character of the spectrum are 
not destroyed when the long range part of the interaction is removed. However, there are 
quantitative changes in the asymptotic behavior of tunneling and vibrational excitations at 
low densities. For the lowest densities (equivalent to r s >, lOOae) the formation of a Wigner 
molecule is observed. In this limit where the long range part of the interaction is crucial, g(x) 
vanishes within finite intervals between the maxima. The structure in the density-density 
correlation function becomes more pronounced with increasing influence of the interaction 
but is less sensitive to total electron spin than expected within Hartree Fock approximation. 
It is worth remarking that the best numerical estimate to our knowledge for the critical value 
of r s for the formation of a Wigner cristal in 3D is also m IOOob |17 |. 

Experiments are frequently performed on 2D quantum dots that are based on GaAs- 
AlGaAs heterostructures. In many of these the electron density corresponds to the inter- 
mediate regime. For the geometry and the electron numbers in typical quantum dots 0, [| 
r s « 3ae (area of the dot w 10 5 nm 2 , number of electrons ~ 10 2 , effective mass ~ 0.07m e , 
dielectric constant « 10). Even for this relatively high electron density in comparison with 
more recent experiments |j| the charge density distribution cannot be expected to be ho- 
mogeneous if we assume that the qualitative features of the electronic properties discussed 
above apply also to the 2D case. In small quantum dots that are weakly coupled to leads and 
that contain only a finite number of electrons, like the experimentally investigated double 
barrier systems ||, a granular charge density should lead to characteristic features in linear 
transport properties as a function of the gate voltage Vq ■ A variation of Vq not only changes 
the charge density or the electron number inside the dot but probably simultaneously its 
shape. In [|] the gate electrode is positioned besides the dot region. For this geometry the 
influence of Vg on the shape of the dot area is obvious. The resonance condition for the 
appearance of a conductance peak in linear transport E-p = En+i — En may be fulfilled more 
than once for given N when Vq is increased. Here Ep is the chemical potential of the 
leads (where no electron-electron interaction is assumed) and En the iV-electron ground 
state energy of the quasi-isolated dot. Such a model could account for the experimentally 
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observed fine structure in the conductance peaks in the linear transport regime reported in 
U which seem to be presently not yet very well understood. 

This work was partly supported by grants of the Deutsche Forschungsgemeinschaft (AP 
47/1-1) and EEC (Science Contract No.SCC*- CT90-0020). 



Figure 1: (a) One-particle density g(x) for N = 3 and (b) N = 4 for different L 
( O.laB < L < 945ae , M = 13). The normalisation is such that / dx g(x) = N . 
When L > lae N peaks begin to emerge. For L >, f 00a B the peaks are well separated. 
Inset: g(x) for a pair potential with an exponential cutoff ( a = lOa^ 1 ) for N = 3 
and L = O.lae (dotted line); L = 9.5a B (dashed-dotted line); L = 472ob (dashed 
line); L = 1417ob (solid line) ( M = 11 ). The minima at large system lengths are 
less pronounced as for a = 0. 

Figure 2: Influence of a pair potential on the energy difference between the lowest 
two multiplets Q (vibrational excitation). The ratio Q(a = 10a I3 1 )/fi(a = 0) is 
plotted versus the particle distance r s for N = 2 (x) and N = 3 (o). (M = 11). 
Inset: logarithm of the energy difference A between the ground state and the first 
excited state within the lowest multiplet versus system length L for iV = 2, a = 
(+) and a = Wa^ 1 (o); iV = 3, a = (x) and a = Wa^ (Y). (M = 11). 

Figure 3: Contour plot of the pair correlation function g c (x,x') (cf eq. [5]) for the 
N = 3 ground state ( S = 1/2 ) at L = 9.5a B (M = 13). 
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